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ABSTRACT 



Aims. Our goal is to study gas properties in large-scale molecular outflows and winds from active galactic nuclei (AGNs) and starburst 
galaxies. 

Methods. We obtained high-resolution (I. "55 x I. "28) observations of HCN, HCO+, HNC 1-0 and HCjN 10-9 of the ultraluminous 
galaxy (ULIRG) Mrk 231 with the IRAM Plateau de Bure Interferometer 

Results. We detect luminous emission from HCN, HCO+ and HNC 1-0 in the QSO ULIRG Mrk 231. All three lines show broad line 
wings - which are particularly prominent for HCN. Velocities are found to be similar (x +750 kms"') to those found for CO 1-0. 
This is the first time bright HCN, HCO^ and HNC emission has been detected in a large-scale galactic outflow. We find that both 
the blue- and red-shifted line wings are spatially extended by at least 0."75 (>700 pc) in a north-south direction. The line wings are 
brighter (relative to the line center intensity) in HCN than in CO 1-0 and line ratios suggest that the molecular outflow consists of 
dense (n > 10* cm"^) and clumpy gas with a high HCN abundance X(HCN)> 10"**. These properties are consistent with the molecular 
gas being compressed and fragmented by shocks in the outflow. Alternatively, HCN is instead pumped by mid-IR continuum, but we 
propose that this effect is not strong for the spatially extended outflowing gas. In addition, we find that the rotation of the main disk, 
in east-west direction, is also evident in the HCN, HCO^ and HNC line emission. An unexpectedly bright HC3N 10-9 line is detected 
inside the central 400 pc of Mrk 231. This HC3N emission may emerge from a shielded, dust-enshrouded region within the inner 
40-50 pc where the gas is heated to high temperatures (200 - 300 K) by the AGN. 



Key words, galaxies: evolution - 
molecules 



1. Introduction 



- galaxies: individual: Mrk 231 — galaxies: active — quasars: general — radio lines: ISM — ISM: 



Mergers and interactions between galaxies trigger massive 
starbursts and feed the growth of supermassive black holes 
(SMBH) in their centers. Powerful outflows have been pro- 
posed to regulate both star-formation and the SMBH growth 
through nega tive feedback i n young galaxi e s (e.g.lHopkins et all 
2009; Murra v et al.1 l2005h . iMurrav et all (l20()5h suggest that 
momentum-driven winds and outflows are the underlying mech- 
anism behind the Faber-Jackson relation for large elliptical 
galaxies. 

In the past years, few studies have targeted outflowing 
molecular gas (e.g. Nakai et al. 1987; Garcia-Burillo et al. 2001; 
IWalter et al.ll2002t ISakamoto et alJI2006i) but recently evidence 
is mounting for massive molecular ou t flows in AGN s and 
starburst gal axies (ISakamoto et al.l 2010t iFerug lio et alj_ 2010t 
IChung et al.1 1201 U lAlatalo et al.1 1201 ll; ISturm et alj l20Tlh and 
high-velocity molecular outflows have been found in OH- 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain) 



megamasers (e.g. lBaanll2007h . It is important to study outflow- 
ing molecular gas because this process removes gas immedi- 
ately available for star formation from the galaxy which impacts 
the evolution of the galaxy on short time scales. Furthermore, 
the molecular gas properties hold important clues to the nature 
and evolution of the outflow. Simulations of molecular cloud 
properties in galactic outflows have been carried out by e.g. 
Narayanan et al. (2008) . 

The ultraluminous infrared galaxy (ULIRG, log(LiR) = 
12.37) Mrk 231 is often referred to as the closest infrared quasar 
(QSO). It is a major merger and hosts powerful AGN activity as 
well as a young, dusty starburst with an extreme star-formation 
rate of ^ 200 M^yr-^ (iTavlor et al.1119991: iGallagher et alJl2002l: 
Lipari et al. 2009). Massive amoun ts of molecular gas re- 
side in an east-west r otating disk dBrvant & Scovilla 119961 : 
iDownes & Solomonll998i) . Radio continuum obser vations show 
jets and outflows ranging f rom pc to kpc-sca le ( Carilli et"al] 
119981: lLonsdaleetai]l2003h . iRupke & Veilleuxl (l2oT jj) discuss 
the various outflows of Mrk 231 and using high resolution 
spectroscopy, they describe a wide, kpc-scale, high-velocity 
1000 km s"') outflow seen in neutral gas (Na I D) absorption. 
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Furthermore, they find that north of the nucleus the radio con- 
tinuum jet is accelerating the outflow to even higher velocities 
1400 kms"') and 3.5 kpc south of the nucleus a slower, 
starburst-driven wind is found. They suggest that the high mass 
outflow rate of Mrk 231 (which is 2.5 times greater than the 
star-formation rate) is a result of the negative feedback from the 
AGN. Interestingly, the outflow has been detected in molecular 
gas manifested i n the broad (75 km s ~') wings in the CO 1-0 
line, as found by iFeruglio et al.l These authors estimate 

the outflow of molecular gas to 700 M© yr ' which coul d empty 
the r eservoir of molecular gas within 10 Myr. Ferugl io et al.l 
(1201 0) argue that the kinetic energy from the supernovae induced 
by the circumnuclear starburst disk is not sufficient to expel 
the gas and that acti on from the AGN is required. Furthermore, 
iMurrav et al.l ( l2005h find that the mass loss rate in a starburst- 
driven outflow does not significantly exceed the star-formation 
ra te (following their equ ation (13)). Th e high mass outflow rates 
of iFeruglio et all (|2010|) and iRupke & Veilleux (.201 1.) seem to 
support the n otion of an AGN-dr iven outflow. This conclusion 
is disputed bv lChung et alj ( 1201 lb . who find high-speed (1000 - 
2000 km s ') molecular outflows in ULIRGs and argue that they 
can be powered by starbursts. They suggest that the mass out- 
flow rate of Mrk 23 1 is overestimated if the CO emission in the 
outflow is optically thin. Clearly, a more detailed study of the 
properties of the gas in the outflow is essential for a better esti- 
mate of the mass in the outflow, but also to understand how the 
molecular gas is carried out of the galaxy and by what process. 
Is the molecular gas entrained in outflowing hot material - or is 
it being expelled by radiation pressure from an AGN and/or the 
momentum flux by supernovae? How can the the gas survive in 
molecular form - and for how long? 

In our IRAM 30m 3 mm EMIR spectrum (ICosta^hola et all 

l20Tlb of Mrk 231 we found that the HCN and HCO"' line widths 
were surprisingly broad, leading us to suspect that the outflow 
signature of Mrk 23 1 could be even more powerful in the rota- 
tional transitions of these high-density tracers. To confirm this 
discovery and to map the extent of the line emission, we asked 
for IRAM Plateau de Bure DDT time to observe HCN, HCO+ 
and HNC in the B-array. 

In sections|2]and[3]the observations and their results - bright 
HCN, HCO^ and HNC 1-0 emission in the outflow line wings 
and a serendipitous detection of the HC3N 10-9 line - are de- 
scribed. In sections [TTI and ITTI we discuss the spatial extent of 
the line wings, and in sections l4n and potential underlying 
causes for the enhanced wing emission are outlined. Outflow gas 
properties and the nature of the outflow are discussed in section 
14.31 In section l4~4l we discuss the origin of the HC3N line emis- 
sion, and finally in section l43] a brief future outlook is outlined. 

2. Observations 

The observations were carried out with the six-element IRAM 
Plateau de Bure array on March 8, 201 1, in B-array. The phase 
center was set to q'=12:56:14.232 and 5= 56:52:25.207 (J2000). 
Mrk 231 was observed for hour angles 2 (i.e. starting 2 hrs 
after source culmination) to 5 and the radio seeing was 0."8. 
The receivers were tuned to a frequency of 85.3125 GHz, the 
center between the two red-shifted line frequencies for HCN 
(y=85.046 GHz) and HCO+ (v^ 85.579 GHz) Zero velocity 
refers to the redshift z=0.042170 dCarilU et al.lll998l) through- 
out the paper, resulting in a spatial scale of 1"=870 pc. We used 
the WideX correlator, which provides a broad frequency range 
of 3.6 GHz, which also included the HNC 1-0 and HC3N 10-9 
lines. 



The bandpass of the individual antennas was derived from 
the bright quasar 3C273. The flux calibration was set on 3C273. 
The close-by quasars 71259 -1-516 (~ 0.32 Jy at 3 mm) and 
71300 + 580 (~ 0.19 Jy) were observed regularly every 25 min- 
utes for complex gain calibration. 

After calibration within the GILDAS reduction package, 
the visibility set was converted into FITS format, and im- 
ported into the GILDAS/MAPPING and AIPS packages for ad- 
ditional imaging. Data were rebinned to 20 MHz (68 km s"' ) and 
cleaned with a robust weighting of 3, resulting in a beam size 
of 1. "55x1. "28 and position angle PA=57°. The corresponding 
rms noise is 0.48 mJy channel"'. 

3. Results 

3.1. HCN 

Spectra of all four detected lines are presented in Fig. [T] 
Prominent line wings are detected in the HCN 1-0 line out to 
velocitie s of 750 km s~', similar to those found in CO 1-0 by 
Feruglio etan(l2010h . If we define the line-center-to-wings ratio 
(LCW) as the ratio between peak intensity at zero velocity to the 
intensity at -1-500 kms 'in the red-shifted line wings, then the 
LCW f or HCN 1-0 is 27 mJ y/3 mJy = 9. In the CO 1-0 spec- 
trum of'Feruglio et al.' (2010) the LCW is 320 mJy/8 mJy = 40. 
Another way of definin g the LCW is to adopt the Gaussian fits in 
Fig. 1 of lFerugfio et al J (12010 ). where two components, one with 
5y(FWHM)= 1 80 km s"' and one with (5y(FWHM)=870 km s"' , 
are fitted to the line center and outflow respectively. Fixing the 
linewidths of the two Gaussian componets to our data (see Fig. [3] 
for the Gaussian fit), we find an LCW for HCN of 7.7 and for 
CO 1-0 we estimate an LCW of 20.6. 

Because the line wings are most prominent for HCN, we 
elected to show the line integrated intensity of the wings for 
this molecule only. The HCN 1-0 integrated intensity maps 
of the line wings and the total line emission are presented in 
Fig. |3] peak- and integrated intensities in Tab. [T] The total line 
emission is beginning to be resolved in the 1. "55x1. "28 beam 
with features extending out to the west and northeast. These are 
ge nerally consistent with sim ilar features in the CO 2-1 map 
by i Downes & SolomonI (fT998.) and in the HCO+ 4-3 map by 
Wil son etalJ (120081) . The central velocity gradient from the main 
disk is evident in HCN - showing the ± 150 kms 'gradient 
seen in CO (Downes & Solomon 1998). The HCN 1-0 velocity 
field of the central component (wings excluded) is presented in 
Fig-El HCN position velocity diagrams along the rotational axis 
(PA=90°) and along the wings (PA=0°) are presented in Fig.|5] 

3.1 .1 . Position of line wings and continuum subtraction 

Because the line wings are broad and relatively faint, the de- 
termination of their positions depends on a reliable continuum 
subtraction. We used two methods of subtracting the continuum 
to test the relative errors in the wing position: subtraction in the 
uv-plane and in the image plane. The continuum was subtracted 
with a zeroth or first order fit to the line-free spectral channels. 
In general we find that the uncertainty in the wing peak position 
is 0."4. The continuum map is presented in Fig. |2l A Gaussian 
fit to the continuum source gives a flux density of 25.0 + 0.6 mJy 
and an upper limit to the source size of 0."39. 

The blue- and red-shifted HCN line-wings appear shifted by 
0."2 - 0."4 (175 - 350 pc) to the southwest of the nucleus as 
defined by the peak position of the integrated HCN 1-0 emis- 
sion (see Fig. [3]). Given the positional uncertainties, we may con- 
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Fig. 1. Plateau de Bure spectra of HCN, HCO^, HNC 1-0 and HC3N 10-9. In the right panels we have zoomed-in on the base 
of the line to show the line wings more clearly. In general, red solid lines show Gaussian fits to the line center (core) line widths 
and are 197 km s 'for HC3N and HNC, 213 km s 'for HCO+ and 236 km s 'for HCN. However, in the top right panel of HCN 
1-0 we instead show two Gaussian fits (one with (5y(FWHM)=180 km s"' and one with 5y(FWHM)=870 km s"'), as discussed in 
section im Zero velocity was set forredshift z=0. 042170. 

elude that the wings have their peak brightness within 0.5 kpc of 
the nucleus. Furthermore, a faint blue-shifted wing feature can 
be seen tentatively 2"to the west of the nucleus (a similar fea- 
ture is seen around the systemic velocity). Faint red-shifted wing 
emission is curving out 2"northeast of the center These features 
should be confirmed at higher sensitivity. 

3.2. HCO^, HNC and HC3N 

In the HCO^ and HNC spectra (see Fig.[T]and Tab.[TJ the wings 
can also be seen - although they are weaker than for HCN. 
The peak (V-0 kms ') line intensity ratio between HCN and 
HCO^ is 1.3, while the ratio in the wings (summed from -1-200 to 
-1-700 km s^'and from -200 to -700 km s"') is ^2.2. For compar- 
ison, the integr ated line ratio is reported to be 1.8 in the single- 
dish survey of ICostagliola etalj (1201 ll) . The peak HCN/HNC 
line intensity ratio is 1.8 and the line wing ratio is a;2.8. For 
comparison, the integ rated line HCN/HNC ra tio is ^2.6 in the 
single-dish survey of ICostagliola et al] (1201 lb . The rotation of 
the main disk is evident for both HCO^ and HNC (see Fig. [T). 



Table 1. Line flux densities". 



Line 


Peak 


Wing* 


Integrated 




[mjy beam"'] 


[mJy beam"'] 


[Jy kms"'] 


HCN 1-0 


28 ± 0.5 


3 


8.4 


HCO+ 1-0 


22 ± 0.5 


2 


5.6 


HNC 1-0 


15 ± 0.5 


1 


3.7 


HC3N 


3 ±0.5 




0.9 



a) The 3mm continuum position is a: 12:56:14.230 ± 0.0013 6: 
56:52:25.206 ± 0.01 (J2000). A Gaussian fit to the continuum source 
gives a flux density of 25.0 ± 0.6 mJy and an upper limit to the source 
size of 0."39. The V LBI radio continuu m peak is: a: 12:56:14.248 6: 
56:52:25.28 (J2000) jCarini et alJl998h . 

b) From the red-shifted wing at V= -1-500 km s"'. 



The line center integrated intensity maps of HCO^ and HNC are 
presented in Fig.|6] 



3 



S. Aalto et al.: Detection of HCN, HCO+ and HNC in the Mrk 23 1 molecular outflow 




Fig. 3. Integrated intensities of the HCN 1-0 line emission. Left panel: The red-shifted wing emission from 350 to 990 kms~'. 
The grayscale range is from to 0.7 Jy kms"' (beam and the contours start at 0.098 Jy kms"' (beam)"' with steps of 0.0998 
Jy kms '(beam) Middle panel: The blue-shifted wing emission from -990 to -350 kms '. Same grayscale and contours as for 
the previous panel. Right panel: The total (from -990 kms 'to 990 kms ') integrated HCN 1-0 line emission from Mrk 231. The 
grayscale range is from to 8 Jy kms"' (beam)"' and the contour levels are 0.45 Jy kms"' (beam)"'x(l,3,5,7,9,ll,13,15,17,19). 
The thick white crosses mark the peak of the integrated HCN 1-0 line emission and the thinner crosses indicate the VLBI position 
(see footnote to Tab.[T]) and 1"=870 pc. 
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Fig. 5. HCN 1-0 position-velocity (pV) diagrams. Left: North-south cut along the minor axis of the disk rotation. The extended 
nature of the line wings is evident here. Center: East-west cut (along the major axis of the rotating disk) and Right: Zoom-in on the 
disk rotation.(l"=870 pc) 



Emission from the HC3N molecule is found for the first time 
in Mrk 231 through the detection of the HC3N 7=10-9 line (see 
Fig. |6]). Its line width is similar (^ 197 kms"') to that of the 
HNC line center. The HCN I-O/HC3N 10-9 peak line ratio is 
7.7 (in the ICostagliola et al.l (1201 lb single-dish survey a lower 
limit of 14 is found). The HC3N 10-9 emission is unresolved in 
the 1. "55x1. "28 beam. 



4. Discussion 

4. 1 . The extent of the HCN line wings 

We confirm the result of lFeruglio et al. l (l20Toh that the line wings 
are emerging from a spatially extended structure. This supports 
the notion that they originate in an extended outflow - rather than 
representing the extremes of a Keplerian nuclear disk. (The no- 
tion of extended Keplerian rotation fails on the extreme required 



enclosed mass.) A two-dimensional Gaussian fit to the integrated 
intensity map of the red-shifted wing results in a full-width-at- 
half-maximum (FWHM) size of 0."85 + 0."2 (l"=870pc). The 
corresponding map for the blue-shifted wing has lower signal- 
to-noise and a Gaussian fits results in a FWHM size of 0."75 + 
0."3. 

Comparing the HCN and HCO^ wings w ith the structure of 
the at omic (neutral) gas as seen in absorption (iRupke & Veilleuj 
1201 a the results are consistent with an almost face-on nuclear 
wind. Rupke & Veilleux (201 1) find the most highly accelerated 
atomic gas in the northern region (with peak velocities occurring 
around 2. "3 (2 kpc) north of the nucleus) where the radio jet is 
interacting with the neutral material. Rup ke & Veilleuxl (1201 lb 
suggest that the jet is entraining material from the large-scale 
wind that accelerates it even more which is a phenomenon seen 
also in other AGNs ([Matsushita et alj|2007h . There is a hint of 
this feature at the 1.5 -2cr level in our HCN integrated inten- 
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Fig. 2. 85 GHz continuum image of Mrk 23 1 , the VLBI position 
is marked with a cross (see footnote to Tab.[TJ. 
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Fig. 4. HCN 1-0 velocity field (including only velocities +200 
kms ') showing the rotation of the molecular disk. The HCN 
1-0 emiss ion rotates in the same east-west fashion as the ^^CO 
emission dBrvant & Scovilldll996l: lOownes & Solomonlll998h . 
The grayscale ranges from -100 to 100 kms~'and the contours 
are 20 kms"' x (-4,-3,-2,-1,0,1,2,3,4). Zero velocity was set to 
redshift z=0.042170 and 1"=870 pc. Cross marks the VLBI ra- 
dio continuum position (see footnote to Tab.[Tli. 



sity map of the red-shifted wind - but this requires confirmation 
at higher sensitivity. Furthermore, the velocities of the nuclear 
wind described by Rupke & Veilleux (201 1) and the HCN wings 
generally agree well. The molecular HCN wings we find occur 
equally close (possibly even closer) to the nucleus as the atomic 



wind. We find no corres pondence to the starburst wind found 
3"south of the nucleus bv'Rupke & Veilleuxl(l201 l\i and we con- 
clude that the HCN wings originate in the nuclear wind, as does 
the atomic wind. 

4.2. Why is the HCN emission enhanced in the line wings? 

The CO 1-0 wing peak intensity is 8 mJy (iFerugho et al ] l2010h 
and for HCN 1-0 it is 3 mJy. For the same source size - and 
transferring to temperature scale - the 7'B(CO)/rB(HCN) 1-0 ra- 
tio in the wings is 1.6 - which is lower than for the line peak at 
center velocity. This is either because the outflowing molecular 
gas is dense, with elevated HCN abundances, or it is caused by 
mid-lR pumping of HCN. 

4.2.1 . Dense molecular gas in the outflow 

If we assume that the HCN molecule is collisionally excited, the 
density in the outflow has to be at least « = lO'' cm for the line 
to begin to be excited, even if the HCN abundance is enhanced . 
We can use a simple non-LTE RADEX dvan der Tak et al ] l2007l) 
model to investigate the CO/HCN 1-0 line ratios in a few possi- 
ble scenarios. The models are not unique and serve only as illus- 
trations of possibilities. (Note also that without multi-transition 
HCN information we are unable to constrain the properties of 
the dense gas in the outflow.) 

We adopt a temperature of t he outflowing gas of Tj^- 100 K 
(the lowest gas temperature fit of iGonzalez- Alfonso et al.l ( l2010h 
from Herschel observations of the molecular gas of Mrk 231) 
and note that for densities n = 10'* cm"^ solutions require an 
HCN abundance (relative to H2) of X(HCN); 10"^ It is actuaUy 
the relative CO and HCN column densities that are constrained 
by RADEX - this means that we have to assume a standard CO 
to H2 abundance ratio (in this case of 10"**) to obtain the HCN 
abundance X(HCN). Such an HCN abundance would be consid- 
ered very high (but see Sect. l4.2.2l below for a brief discussion). 
For a number density « = 10^ cm"^ a line ratio between CO and 
HCN 1-0 of 2 can be achieved for optically thin (r < 1) CO and 
HCN with a column density ratio (A^(HCN)/A^(CO)) of 5 x 
If we assume a CO abundance of 10"^ (relative to H2), this would 
imply an HCN abundance of X(HCN)=5 x 10 **, which is high 
- but not unusually so (see Sect. 14. 2."2l below). Intrinsic bright- 
ness temperatures for optically thin emission would be low, how- 
ever, and the emission would have to have a large filling fac- 
to r The low HCN/H'^CN 1-0 intensity ratio of 5.5 observed 
bv lCostagliola et aTl (|201 ih suggests that the HCN 1-0 emission 
towards Mrk 23 1 is not optically thin. This is a global ratio, how- 
ever, and the ratio for the wing component is difficult to deter- 
mine. 

A low intensity ratio between CO and HCN can also be 
achieved by a clumpy outflow where the individual clumps are 
dense and optically thick. A rB(CO)/rB(HCN) ratio of ^2 can 
be obtained for a range of abundances depending on how high 
the CO column density can become in each clump. If we avoid 
the most extreme optical depths in CO, then HCN abundances of 
X(HCN)= 10 - 10"^ win fit the observed Une ratio. This sce- 
nario (or the optically thin one above) does not allow any CO 1-0 
emission to emerge from a diffuse, low-density (n = 10^ - lO"' 
cm ■') component in the outflow - such a component would only 
contribute to the CO 1-0 luminosity, not the HCN. Therefore, 
if the molecules are collisionally excited, the low intensity ra- 
tio between CO and HCN in the outflow leaves little room for a 
low-density molecular component. 
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4.2.2. Outflow abundances 

HCN abundances: Elevated HCN abundances have been ob- 
served in outflows in the Galaxy. For example, HCN abundances 
of 14 X 10"*^ are found in the outflow lobes of region s of mas- 
sive star-formation (e.g. lSu et al BoOTI : iLiu et al The fact 
that the outflow in Mrk 23 1 is seen (tentatively) in H2O absorp- 
tion dGonzalez- Alfonso et alJl20Tol) may suggest that shocks are 
involved in removing the water from the grains. Are we seeing 
indications of HCN abundances being enhanced in the shocks of 
the outflow? Ev i dence that this can be the case are presented 
by iTafaUa et alJ (l20Toh . They argue that HCN is mainly pro- 
duced in shocks by the reaction of N with CH2 - where CH2 is 
formed in the reaction between C and H3 or C^ with H2. Thus, 
for HCN to be selectively produced in the shock, a sufficient 
abundance of free C is required. SiO emission from shocks in 
galaxy-scale outflows has been observed in the starburst wind of 
M 82 (jGarcia-Burillo et al. 2001) and may also be expected to 
be elevated in the Mrk 23 1 outflow. 

Note also that a larg e HCN abundance is a general re- 
sult of high temperatures. iHarada et al.l (1201 Ol) have shown that 
HCN abundances may reach 10"^ in regions where r!300 K. 
They also find enhancements of H2O abundances in their high- 
temperature models - therefore HCN and H2O abundance en- 
hancements may be a result of high temperatures and not specif- 
ically shock chemistry. A potential scenario is that the molecular 
abundances reflect the high temperatures near the AGN and re- 
tain these abundances in the outflow. However, without a source 
of heating, the outflowing gas will cool very quickly if it re- 
mains at high density and potential ion-neutral reactions (involv- 
ing e.g H^, Hj , O^, HsO^ etc) may destroy HCN on short time- 
scales. So it is unclear whether the gas can retain relic high- 
temperature abundances of HCN. Furthermore, at least part of 
the outflow will be exposed to the X-rays of the nuclear AGN, 
which means that it could possibly influence the chemistry of the 
outflowing gas (e.g. Malonev et al. 1996; Meiierink & Spaans 
QoOStlMeiierink et alJl2006l) . ' 

HNC and HCO^ abundances: The HCN intensity is greater 
than that of HCO^ by a factor of two in the line wings. This 
could be caused by excitation or by HCN being more abun- 
dant than HCO+. In Galactic outflow sources, I(HCN)>I(HCO+) 

i and abundance rat ios HCN/HCO^= 10-50) are not uncommon 
Tafafla et a DIMS and may simply be a result of the enhanced 
HCN abundance in the outflow due to shocks and high temper- 
atures. The line wing HCN/HNC ratio of 2.8 is within th e range 
of wh at is typically observed in galaxies (e.g. Costagliola et alj 
1201 Ih . In general HNC is exp ected to be destroye d in shocks - 
and at high temperatures (e.g. lSchiIke et alJll992h . For the sec- 
ond model in section 14.21 an abundance ratio HCN/HNC and 
HCN/HCO^ of a!20 may reproduce the observed 1-0 line ra- 
tios in the line wings. Both these abundance ratios are consistent 
with values found for Galactic shock regions. 



4.2.3. IR pumping 

An alternative scenario to collisional excitation of HCN is 
pumping by 14 yum mid-IR continuum. The rotational levels 
are then populated via the IR bending modes. Vibrational line 
emission from HCN has be en observed in ULIRGs such as 
Arp 220 dSalter et alJ I2OO8I) and in LIRGs like NGC 4418 
(ISakamoto et alJ 20101) . showing that the HCN molecule is 
pumped by mid-IR absorption in these galaxies. For HCN 
to be pumped a minimum 1 4 /zm brightness temperature of 
rB(IR)=85 K is required (e.g. lAalto et"d]l2007h . If IR pumping 



dominates the HCN excitation no dense clumps in the outflow 
are needed to explain the HCN emission. It is therefore impor- 
tant to settle the question of the excitation of HCN in the outflow. 
Furtherm ore, both HCO^ a nd HNC may be pumped by mid-IR 
emission jAalto et al.ll2007h through bending modes in the same 
way as for HCN. HNC is pumped already at 21 and HCO^ 
at 12 jim. 

In the SED for Mrk 231 (iGonzalez-Alfonso et al.ll2010l) the 
14 jum regime is dominated by a hot component with Td^lSO- 
400 K emerging from the inner 40 - 50 pc. One would expect the 
strongest pumping of HCN close to this hot component, i.e., the 
nucleus. Furthermore, the plane of the disk is likely opaque at 
14 jum, but the gas lifted out of the plane in the outflow will have 
an unobstructed view of the nuclear region. Very tentatively, one 
might argue therefore that the HCN in the wings is pumped more 
strongly than that in the core of the line. However, the extended 
nature of the wings implies that the gas is at least a few hundred 
pc from the IR source and the resulting brightness temperature 
produced by this source is going to be well below 85 K, render- 
ing the pumping scenario unlikely. This argument holds also for 
HNC and HCO+. 

It is interesting to compare the potential HCN pump- 
ing scenario to the case fo r water detected by Herschel 
dGonzalez- Alfonso et alJl20To[) . H2O is pumped at significantly 
longer wavelengths (than HCN) from 45 to 75 and may 
therefore be pumped by the warm, mo re extended dust compo- 
nent (r=120 pc disk with Td = 95 K dGonzalez-Alfonso et al.l 
"201^). This region is smaller than that of the HCN line wings 
and the H2O absorption line profile also appears different from 
that of the HCN line wings. It is therefore possible that H2O is 
IR-pumped while HCN is not. 

4.3. Gas properties and the nature of the outflow 

Our interpretation of the HCN-luminous wings of Mrk 231 is 
that the molecular gas in the outflow is clumpy and dense, with 
an enhanced HCN abundance. The similar spectral shapes (and 
spatial extent within the errors) of the CO and HCN 1-0 line 
wings suggest that the emission is emerging from the same out- 
flowing gas. This would then imply that there is no significant 
evolution (e.g. star-formation or evaporation) of the cloud prop- 
erties along the molecular part of the outflow. This is inconsistent 
with the simulations of molecular cloud properties entrained in 
outflows by Narayanan et al, (2008), who predict a shorter life- 
time and smaller spatial extent of the dense gas phase - as well 
as ongoing star-formation in the outflowing molecular gas. The 
Mrk 23 1 molecular outflow appears instead to be dominated by 
the dense phase without much low-density molecular gas. Is this 
because radiation dissociates the low-density gas, or because it 
becomes compressed in the flow? Alternatively, the molecular 
cloud structure largely reflects the interstellar medium structure 
from where it becomes launched into the outflow. Note that the 
spatial extent of the CO and HCN wings must be better deter- 
mined for a deeper exploration of the evolution of the gas in the 
outflow. 

In an AGN-driven outflow, the molecular gas can be en- 
trained b^jijjienm^ wind or driven by momentum depo- 
sition. iRupke & Veilleuxl d201 Ih find that the neutral gas has 
higher velocities than the ionized gas in Mrk 231, which is not 
consistent with entrainment by a hot wind. They sugges t the 
AGN-driven two-phase model by iHopkins & ElvisI d201( ft as a 
possible scenario for the neutral wind of Mrk 23 1 : Close to the 
center, the AGN acts on the hot interstellar medium, powering 
an outflow that causes the cold molecular clouds to expand. 



6 



S. Aalto et al.: Detection of HCN, HCO+ and HNC in the Mrk 23 1 molecular outflow 



T 




Fig. 6. Integrated intensities (from left to right) of the HCO^, HNC 1-0 and HC3N 10-9 line emission. The grayscale 
range is from to 5.3, to 3.3, and to 1.0 Jy kms"' (beam)"' respectively. The contour levels are 0.26, 0.16, and 0.04 
Jy kms"' (beam)-ix(l,3,5,7,9,ll, 13,15,17,19). (1"=870 pc and and zero velocity was set to redshift z=0.042170 for all lines). 
Cross refers to the center of the 85 GHz continuum. The maps are integrated over +200 km s 'and exclude the line wings. 



which makes them more susceptible to being carried outward 
by radiation pressure (and to become ionized). However, this re- 
quires a mechanism to rapidly compress the molecular gas in 
the subsequent radiation pressure-driven wind for it to be con- 
sistent with our observations. Alternatively, radiation pressure- 
driven outflows could be launched directly, close to the black 
hole (e.g. King 2003). Dense molecular gas from the inner re- 
gions of Mrk 23 1 may then be compressed even more in the out- 
flow without being initially disrupted. High HCN abundances 
may indicate shock enrichment in the outflow and/or that the 
chemistry is influenced by the AGN. It is interesting to speculate 
if HCN-luminous wings may be a general signature of AGN- 
driven outflows in luminous galaxies. High-sensitivity HCN sur- 
veys of starburst and AGN-driven outflows will provide us with 
an answer 

We note that Chung et alj (1201 lb detected CN emission only 
toward the AGN-driven objects, in their stacked survey of out- 
flows in ULIRGs, not toward the starbursts. This is either be- 
cause of higher gas densities in these outflows, or because of CN 
abundance enhancements, and also emphasizes the importance 
of studying the molecular gas properties of outflowing AGN and 
starburst galaxies. 

4.4. Serendipitous HC3N detection: warm, sliielded gas in 
fhe nucleus. 

The fitted upper size limit of the HC3N 10-9 integrated line 
emission is FWHM=0."48 (400 pc) and it is centered on the 
nucleus of Mrk 231. Emission is not detected in the line wings. 
The peak-to-peak line ratio between HC3N 10-9 and HCN 1-0 
is 0.14+0.015 - which makes the HC3N line emission (relative 
to HCN) brighter than in, for example, NGC 1068 or NGC 777 1 
- and almo st as bright as in Arp 22 dCostagliola et al.ll201 ll) . 
Models by iHarada & Herbstj (l2010h indicate that HC3N abun- 
dances may be enhanced substantially (orders of magnitude) 
in warm (> 200 K) regions in the midplane of dense disks 
around AGNs. These regions do not qualify as XDRs - X- 
ray dominated regions - but may be buried within the XDR. 
According to models bas ed on the recent Herschel observations 
(IvanderWerf etal.ll2010t) . the gas in the 1 . 1 kpc starburst disk is 
exposed to intense UV emission - which should effectively de- 
stroy the HC3N molecules. Since the size of the HC3N emitting 
region is smaller than the suggested dimensions of the starburst 



disk, the HC3N emission may emerge from a region between the 
inner edge of the starbur st disk and the AGN dom inated XDR. In 
the dust SED model bv iGonzalez-Alfonso et al.l (1201 0) the nu- 
clear hot component with rd=150-400 K originates in the inner 
40 - 50 pc. This may therefore be the location from where the 
HC3N is emerging - consistent with the expected HC3N abun- 
dance enhancement in warm, dusty regions. The high bulk tem- 
perature required for significant abundance enhancement is con- 
sistent with the AGN as the heating agent. It is possible that this 
dusty, warm regio n represents the fi nal stage of the obscured, X- 
ray absorbed (e.g. lPage et al.ll2004l) accretion phase of the QSO. 

4.5. Future outiook 

We have seen that the molecular outflow in Mrk 23 1 has surpris- 
ing gas properties and physical conditions that contain important 
clues to its origin, physical conditions, and chemical composi- 
tion. To fully understand the enhancement of HCN (and also the 
HCO^ and HNC emission), multi-transition, high-resolution ob- 
servations are necessary. Other species such as multi-line CO 
observations, isotopic variants of CO, HCN, HCO^ and HNC, 
shock tracers like SiO and CH3OH will provide essential, new 
pieces to this puzzle. These studies should be accompanied by 
radiative transport and astrochemical modeling. Deep imaging 
of CO and HCN will test the notion that the emission is co- 
spatial and originates in the same gas. 

The detection of HC3N emission toward the inner region of 
Mrk 23 1 raises interesting questions about regions of shielded 
gas near the nucleus of this ULIRG QSO. High-resolution ob- 
servations will pinpoint the exact location of the HC3N emitting 
regions, allowing us to better understand the underlying mecha- 
nisms and evolution of this powerful galaxy. 

5. Conclusions 

In this paper we present the first high-resolution observations 
of HCN, HCO+ and HNC 1-0 emission from the wings of the 
prominent molecular outflow of the QSO ULIRG Mrk 231. The 
velocities (k 750 kms~') are similar to those akeady found for 
the CO-outflow. The line wings are spatially extended by at least 
0."75 (>700 pc) in a north-south direction. This is the first time 
prominent HCN and HCO^ emission has been detected in a 
large-scale galactic outflow. The east-west rotation of the main 
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ratio is 1 .8 and the line wing ratio is 2.8. This is within the range 
of what is typically observed in galaxies. 

An unexpectedly bright HC3N 10-9 line was detected in the 
central 400 pc of Mrk 23 1 . This molecule is often associated 
with young star-forming regions and tends to be destroyed by in- 
tense UV and particle radiation. We find it unlikely that the emis- 
sion originates in the starburst region of Mrk 23 1 and instead 
suggest that the HC3N emission is emerging from a shielded, 
dust-enshrouded region within the inner 40-50 pc where the gas 
is heated to high temperatures (200 - 300 K) by the AGN. We 
speculate that this dusty, warm region represens the final stage 
of the obscured. X-ray absorbed accretion phase of the QSO in 
Mrk 231. 
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Fig. 7. Central HCO+ (top) and HNC (bottom) 1-0 velocity 
fields. The grayscale range is from -64 to 9 km s'^top) and from 
-43 to 24 km s~Hbottom). The contour levels are 5 kms~^x(- 
4,-3,-2,-1,0,1,2,3,4). (1"=870 pc and zero velocity was set to 
redshift z=0. 042170 for all lines). Cross refers to the center of 
the 85 GHz continuum. The emission is integrated over ±200 
km s~^and excludes the line wings. 



disk is also evident in the HCN, HCO"^ and HNC 1-0 narrow 
component line emission. 

Compared to CO 1-0, HCN 1-0 appears relatively enhanced 
in the line wings by factors of 2-5. This suggests that a large 
fraction of the gas in the outflow is dense n > 10^ cm"', which 
is consistent with the molecular gas being compressed and frag- 
mented by shocks. The HCN abundances may be significantly 
elevated (Z(HCN)> 10"^) in the outflow, which can be caused 
by shocks and/or high temperature chemistry. 

Alternatively, the excitation of HCN could be affected by 
mid-IR radiation, but we suggest that this is less likely because 
of the extended nature of the outflow. 

The HCN/HCO"^ Une peak ratio is 1.3 for the line center and 
2.2 in the line wings. l(HCN)>I(HCO^) is not uncommon in out- 
flowing gas in Galactic sources. The peak center HCN/HNC Une 
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